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Abstract: We present a general method to calculate the connected correlation func- 
tion of random Ising chains at zero temperature. This quantity is shown to relate to the 
surviving probability of some one-dimensional, adsorbing random walker on a finite in- 
tervall, the size of which is controlled by the strength of the randomness. For different 
random field and random bond distributions the correlation length is exactly calculated. 
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1 Introduction: 



The perfect crystal is the result of physical abstraction in real materials one should 
always reckon on the existence of impurities and different types of lattice defects. There 
are important problems in the field of random systems such as localization in a disordered 
medium, the spin glass behaviour, diluted magnets etc. In order to obtain a theoretical 
understanding of these phenomena magnetic models with quenched disorder have been in- 
troduced and studied by different methods (for recent reviews see Refs.[l-3]). The random 
magnetic models have unusual low-temperature properties. These are a consequence of 
the complex structure of low energy metastable states generated by quenched disorder and 
frustration, which are beleived the main ingredients of spin-glass behaviour [4]. The sim- 
plest system in this field is the random Ising model, which can be experimentally realised 
as diluted antiferromagnets[5]. The model already in one dimension shows interesting 
features, although its physical quantities are singular only at zero temperature. 

The one-dimensional random Ising model - inspite of its low dimension - is exactly 
soluble only for a few specific form of randomness. One of those is the random bond Ising 
chain in a uniform field (RBIM) defined by the Hamiltonian: 



Here Oi = ±1 and the exchange integral Ji^+i is equal to J > 0, with probability p and 
— J with probability q = 1 — p, and the bond disorder is quenched. For this model the 
ground state energy, the zero point entropy and magnetization have been calculated by 
Derrida and co-workers[6] (see also in Ref[7]). 

Another type of random system is the ferromagnetic Ising model in a random field 
(RFIM) with the Hamiltonian: 



At T = the thermodynamic properties of this model have been calculated for the binary 
distribution of the random fields [6]: 




(1) 




(2) 



P(h) = p5(h -H)+ qd(h + H) 



(3) 



as well as for the asymmetric binary distribution[8]: 




(4) 
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The first exact results at T ^ have been obtained by Grinstein and Mukamel[9] with the 
diluted symmetric binary distribution: 

P(h) = Z[6(h -H) + S(h + H)] + qS(h) (5) 

in the limit H — > oo. Very recently the non-linear and higher order susceptibilities of the 
model have also been determined [10]. 

A series of RFIM-s with continuous random field distribution has been studied by Luck 
and Nieuwenhuizen[ll-13] at arbitrary temperatures. One of these models is characterised 
by the diluted symmetric exponential distribution R(x)dx as: 

R( x ) = l e ~ x + qS{x) (6) 

with hi = Hxi, H > and — oo < Xi < oo. 

As far as the correlation functions of random Ising models are concerned only a few 
results are available. Derrida[14] has pointed out that the spin-spin correlation function is 
not a self-averaging quantity, therefore its average differs from the most probable value[15]. 
Exact results in the whole temperature range are available for the Grinstein-Mukamel 
model[9] as well as for the diluted symmetric exponential distribution in eq(6)[12]. Most 
recently Farhi and Gutmann[16] have calculated the zero temperature correlation function 
of the RFIM with the binary distribution in eq(3). We also mention a related exact study of 
the pair-correlation function of a one-dimensional lattice gas model in a random potential 
at zero temperature[17]. 

In this paper we study the connected correlation function of random Ising chains 
defined as: 

X(l) = [< <Ti<T i+ i > - < <Ti >< a l+ i >] av (7) 

where < ... > denotes the thermodynamic average and stands for the quenched 

average over the random variables. The thermodynamic average in eq(7) could be non- 
zero even at T = provided the ground state of the system is highly degenerate. This 
can be seen in two-dimensional frustrated models without randomness, in which the zero 
temperature correlations either decay as a power low[18-21] or exponentially [22]. 

For random Ising chains the correlation function in eq(7) is usually calculated in 
the transfer matrix formalism [9, 12]. At T = 0, however, one may use another approach 
based on an analysis of the degenerate ground state configurations. For a given quenched 
disorder, due to frustration, there are spins in the system which are "loose", i.e. they are 
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free to point in any direction. These loose spins may form a domain, and the connected 
correlation function is non-zero, only if both spins considered belong to the same domain 
of loose spins. At this point calculation of the correlation function at T = is essentially 
reduced to an investigation of the size distribution of domains of loose spins, which in turn 
is equivalent to a one-dimensional random walker problem on a finite intervall. 

The setup of the paper is the following. In Section 2 we present our method to 
calculate the connected correlation function for random Ising chains at T = 0. In Section 
3 and 4 correlation lengths are calculated for RFIM and RBIM, respectively. Finally, the 
results are discussed in Section 5. 

2 Connected correlations and their relation with random walkers 

In this Section we develope a formalism to calculate the connected correlation function 
for random field Ising models at T = 0. It will be shown in Section 4 how these results 
can be applied for the random bond problem. 

The random fields we consider have a discrete distribution, furthermore the possible 
values of hi are integer multiples of a unit, denoted by H, i.e. hi = m x H. For the 
binary distribution in eq(3) m = ±1, while for the diluted symmetric distribution in eq(5) 
m = 0, ±1. Continuous distributions, like in eq(6) can also be discretised, an example is 
shown in Section 3.3. In this way one can also treat the asymmetric binary distribution in 
eq(4) provided the ratio of the parameters Hi and H is rational. 

The structure of ground state configurations of a RFIM with discrete randomness is 
thoroughly analysed in the literature (see c.f. in Refs[2,3]). In the weak-coupling limit, 
when 2 J < min{\hi\} = h m i n the spins are frozen to the direction of the local fields. For 
stronger couplings, so that 2J > h m i n there is a tendency for neighbouring spins to align 
parallel with each other, so that domains of parallel spins are formed. With increasing 
value of the coupling the average size of a domain increases, but the ground state never 
consists of one single domain. It can be understood, since the necessary energy to create 
a domain wall - 2 J - can be accumulated from fluctuations of the random field, even if the 
local field is arbitrarily small. 

The size of random field fluctuations is characterised by the integrated random field 
function, defined as: H(k) = Yli=i hi- For an illustration we draw this function on Fig 1 
for a given random field distribution together with the corresponding ground state config- 
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uration of the system at some value of the coupling J. As is seen on this Figure the first 
domain wall is located between spins 3 and 4, at a local maximum of H(k). Indeed, the 
sum of random field energies H(6) — H(3) < —2J covers the cost of creation a domain 
wall. 

The position of the second domain wall, however is not unique: it can be at any of the 
three degenerate local minima located at (6, 7), 9, 10) or (11, 12). Energetically there is no 
difference between these configurations, since the corresponding random field energies are 
the same: H(6) = H(9) = H(ll). As a consequence in the intervall {7, 11} the position 
of the spins in the ground state is not fixed. Such a region will be called as a domain of 
"loose" spins (DLS). 

Based on this example we can easily postulate the properties of a DLS. First let us 
restrict ourselves to a DLS which separates a j ... j and a f ... f domain, as shown in our 
example on Fig 1. Such a DLS is bounded by two degenerate local minima. Inside the DLS 
the integrated random field function relative to its value at the boundaries - AH(k) - does 
not exceed 2 J, thus the energy necessary to create a domain wall can not be accumulated 
from the random field. On the other hand in both directions outside of the DLS AH(k) 
exceeds 2 J, before crossing zero. This last condition ensures the existence of j ... } and 
| ... | ferromagnetic domains at two sides of the DLS. 

The other type of DLS separating | ... f and | ... | domains can be characterised 
similarly. In this case after the transformation Oi — > — Uj, hi — > —hi the previous consid- 
erations can be applied for —H(k). We mention, if 2J/H is an integer the degeneracy of 
the ground state is higher, than for a slightly larger or smaller value of the coupling. It is 
connected to the fact that in this case conditions both for f| and || DLS can be satisfied 
at the same time. In the following we are not going to deal with such situations, thus our 
considerations apply for 2J/H ^ integer. 

Since position of spins in a DLS is not fixed these regions are the source of non-zero 
ground state entropy in a RFIM. These regions are also responsible for non-vanishing value 
of the connected correlation functions at T = 0. It is easy to see that the thermal average 
of x(l) m e q(7), which is now performed over the degenerate ground state configurations, 
is no n- zero only if the two spins are in the same DLS. If one spin is in a ferromagnetic 
domain its value is the same for all ground state configurations, consequently the con- 
nected correlations are zero. On the other hand if two spins are in different DLS they are 
independent variables, thus again the connected correlations are vanishing. 
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The above relationship between x(Z) and DLS makes it possible to calculate the leading 
behaviour of the connected correlation function in a RFIM in a simple way. To do this 
one should i) first consider DLS regions of length / > / and determine the W(l) probability 
that a point of the line belongs to one of those DLS. ii) Then find the probability that the 
other reference point of x(Z) is also on the same DLS and iii) finally calculate the thermal 
and quenched averages in eq(7) with the condition that both endpoints of x(Z) are on the 
same DLS. In the following we show that the leading behaviour of x(Z) is determined by 
the probability W(l), which has an exponential dependence on Z, whereas the probabilities 
indicated in ii) and iii) are comperatively negligable, since they depend on Z as power laws. 

To calculate the probability W(l) one can use a geometrical interpretation of a DLS 
as a one-dimensional random walker with steps hi/H on an intervall consisting of 

h\ +1 (8) 

points. Here [x] denotes the integer part of x and x is non-integer. The walker starts at one 
endpoint of the intervall and after Z steps made on the strip returns to the same endpoint. 
For large Z the leading Z dependence of W(l) is exponential, and W(l) corresponds to the 
surviving probability of the adsorbing random walker: 

W(l) ~ exp [-Z~/e(L)j (9) 

Next we show that the probabilities mentioned in ii) and iii) have weaker Z dependence 
than the surviving probability. It is clear from a simple geometrical consideration that 
the probability in ii) is at most ~ 1/7. On the other hand to estimate the conditional 
probability in iii) one should first notice that the different ground states are characterised 
by one parameter, the position of the domain wall. Thus an average over this parameter is 
equivalent to the thermal end quenched averages in eq(7). Using the fact that the number 
of possible positions of a domain wall in a DLS of length Z is proportional to Z one can 
estimate the conditional probability in iii) as ~ Z -2 . 

The connected correlation function x(Z) is then obtained by summing the probabilities 
for Z > Z with the result in leading order: 

X il) ~ W{1) ~ exp [-Z/£(L)] (10) 

where L is given in eq(8). In the following Section the correlation length £(L) will be 
calculated for different random field distributions. 
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3 Random field Ising models 



The surviving probability in eq(10) can be most easily calculated in the transfer matrix 
formalism [23]. In this case the elements of the transfer matrix T(n, m), n, m = 1, 2, L 
are given as the probability of a step from position m to n. In the RFIM lanquage 
T(n,m) = P(h(n, m)), where h(n,m) = (n — m) x H. A matrix-element of T is zero, 
whenever the corresponding /i(n, m) is not contained in the set of random fields of the 
model. The leading eigenvalue of the transfer matrix - A(L) - is connected to the surviving 
probability as: 

W{l)~\{L) 1 (11) 
thus the correlation length in eq(10) is given by 

For a general RFIM simple analytical results can be obtained for 2J/H < 1 and in the 
strong coupling limit 2J/H ^> 1. In the former case L = 1, the intervall of the walker 
consists of one single point: 

A(l) = P(0) 2J/H < 1 (13) 

Therefore non- vanishing correlations can only be present in diluted models (see cf. eqs.(5) 
and (6)). In the strong coupling limit, which corresponds to L > 1 we consider distri- 
butions with zero average < hi >= 0. Then the finite size corrections to the leading 
eigenvalue are quadratic: 

1 - A(L) ~ L~ 2 ~ (H/J) 2 H/J < 1 (14) 

which follows from the Gaussian nature of the free random walk [23]. According to eqs.(12) 
and (14) 

£(J/H) ~ ( J/Hf J/H » 1 (15) 

For intermediate (non-integer) values of 2 J/ H the leading eigenvalue of the transfer matrix 
can be calculated numerically, so that - together with the asymptotic relation in eqs(15) 
- one can in principle obtain the correlation length of connected correlations for all types 
of RFIM-s. In the following we present three examples, in which the calculation can be 
performed analytically for all non-integer values of 2J/H. 
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3.1 Binary distribution 



The transfer matrix corresponding to the distribution in eq(3) is tridiagonal and given 



as: 



/0 p 
q p 
q p 



V 



\ 



p 

q 0/ 



The (non-normalized) leading right eigenvector of this matrix is 
and the corresponding leading eigenvalue 



Ai(L) = 2 v ^cos (j-^j 



(16) 



(17) 



(18) 



The correlation length is then can be obtained from eqs(12) and (8). This result agrees 
with that of Ref[16]. 

3.2 Diluted symmetric binary distribution 



The transfer matrix corresponding to the distribution in eq(5) is symmetric and tridi 
agonal: 

( q p/2 \_ 

p/2 q p/2 

p/2 q p/2 



\ 



p/2 
p/2 q / 



(19) 



The leading eigenvector of T 2 is the same as in eq(17), however with q/p = 1 and the 
leading eigenvalue is given as: 



X 2 (L) = q + pcos 



TV 



L + l 



(20) 



For L = 1, i.e. 2J/H < 1, A 2 (l) = q and the general relation in eq(13) is recovered, which 
agrees also with the result of Grinstein and Mukamel[9] as 2J/H — > 0. 



3.3 Diluted symmetric exponential distribution 



We consider the discretised version of the distribution in eq(6) when the allowed values 
of the random field are hi = H x i, i = 0, ±1, ±2, ... and the probability distribution is 
given as: 





h 






H 



5(h - Hi) + c 5(h) 



(21) 



Here a > and the distribution is normalised with c = p/2(expa — 1), Co = q — c and 
p + q = 1. Taking the limit 



a^0+, #->0, # = #/(exp a - 1) = finite 



(22) 



one arrives to a continuum description in the variables hi = Hxi, — oo < X{ < oo, with the 
probability distribution R(x)dx given in eq(6). 

We write the transfer matrix of the problem in terms of the variables u = exp(— a), 
ujq = q and k = c/q as: 



( 1 

KUJ 
KU) 2 



KU KUJ 
1 KUJ 
KUJ 1 



(23) 



\^ L " 1 ... 1 / 

We mention that the same transfer matrix belongs to a directed polymer [23] on a strip 
of width L on the square lattice. In this case ujq and ui are the monomer fugacities for 
steps along and perpendicular to the strip, respectively and k denotes the statistical weight 
corresponding to a turn of the chain by 90°. 

The eigenvalue problem of T 3 in eq(23) is similar to that of the unrestricted directed 
polymer [24]. The leading eigenvector is given as 

L + l~ 



$ 3 (/) = cos 
where <j> is the smallest root of the equation: 

(£-1)0 



l- 



(24) 



cot 



sin i 



COS — UJ 



(25) 



The leading eigenvalue is then: 



> / t \ (1 — uj 2 )kujo , 

A 3 L = >— ^ + 1 - k)uj 

l — 2u cos 4> + uj 2 

9 



(26) 



The correlation length can be obtained from eq(12) using the correspondences between 
u>, cuo, k and the original parameters of the distribution in eq(21). 

Now we evaluate our results in eqs(25) and (26) in the continuum limit of eq(22), 
which reads as w — > 1~. The smallest root of eq(25) is proportional to 1 — uj, so that 
written it in the form of <p = 7(1 — <^) one gets the following relation from eq(25): 

tan - 

Then the leading eigenvalue is given by 

A 3 

in the continuum limit. These results are identical to those obtained by Luck and Nieuwen- 
huizen[12] using the continuous distribution in eq(6) and taking the non-trivial limit 
T — > + . We note that the correlation function calculated with the continuous distribution 
is discontinuous at T = 0, thus the limits H — > and T — > can not be interchanged. 



J_ 



1 



1 + 7 2 



(27) 



(28) 



4 Random bond Ising model 

The formalism developed in Section 2 can also be applied to calculate x(Z) for the 
RBIM with the Hamiltonian in eq(l). It is easy to see that after the transformation Oi — > 
Gihjhi in eq(2) one arrives to a RBIM with J^i+i = J /hihi + \. The inverse transformation 
is given by 




thus the random fields can be expressed through the random bonds as 

^ = *n(^) (30) 

The random fields in eq(30) can take the values ±h like to the binary distribution in eq(3), 
however these h{-s are correlated in different sites, since 

<h l h l+n >=(p-q) n (31) 

Thus the symmetric distribution with p = q = 1/2 is exceptional, in which case the RBIM 
is equivalent to a RFIM with the symmetric binary distribution in eq(3). 
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For p 7^ q the probability is connected with bond variables, thus with the sign of the 
product of two consecutive random fields: 

p hi-xhi > 



Pihi-xhi) = 



q K-xK < 



(32) 



To write down the transfer matrix of this problem we work with the bond variable H (i) = 
(H(i - 1) + H(i))/2, which may take L - 1 different values: 1/2,3/2, ...,L - 1/2. The 
one-step transfer matrix Trb{i, i + 1) is different if the number of the step i is even or odd. 
Therefore one should consider the two-step transfer matrix defined as: 

\ 



qp 



T RB (i,i + 2) 



qp p 
q 2 qp 
qp q 2 



p qp q 



qp p 

2 



qp 
qp q 2 
p 2 qp 



qp 
Q 2 



p* 
qp 



\ 



(33) 



/ 



The transfer matrix connecting the states in steps i + 1 and i + 3 is the transpose of 
TRB(i : i + 2) : thus both have the same eigenvalue spectrum. Writing the leading eigenvalue 
as \ 2 RB the correlation length is then obtained from the logarithm of \rb through eq(12). 
Since the transfer matrix in eq(33) is non-symmetric and non-tridiagonal we could solve 
its eigenvalue problem analytically only in a few special cases. 

In the limit p — > the transfer matrix is tridiagonal in linear order of p and can be 
solved by the same eigenvector as T 2 in eq(19). The leading eigenvalue is given as 



Ai tB = l-|(l-cos^)+0(p 2 ) 
In the symmetric distribution p = q = 1/2 the right eigenvector is given by 



$ RB (2l) = ®rb{21 + 1) = sin 



(21 + 1) 



n 



L + l 



(34) 



(35) 



for I = 0, 1, L/2 — 1. The corresponding eigenvalue 



^rb = cos 



7T 



L + l 



(p = q= 1/2) 



(36) 



is the same as for the RFIM with symmetric distribution in eq(18), which is in accord with 
our previous claim about the equivalence of the two problems. 
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Finally, we consider the limit q — > 0. Then the eigenvector is given in leading order 

as: 

tf(2Z) = *(L-2/ + l) = g^r 

, , (37) 
*(2Z - 1) = *(L - 2/ + 2) = Q 1 " — 

for Z = 1,2, (L — l)/2. The corresponding eigenvalue: 

Xrb = q 1/L q « 1 (38) 

Analysing the p-dependence of the correlation length one can see that it starts with zero 
in the ferromagnetic limit p — > 1, stays finite for non-zero concentration of ferromagnetic 
bonds and finally diverges in the antiferromagnetic limit p — > 0. 



5 Discussion 

In this paper the connected correlation function of random Ising chains is investigated 
at zero temperature. Our study is based on an analysis of the ground state configurations of 
these systems. Due to frustration there are spins in the ground state which are "loose", i.e. 
their position is not fixed by the interaction and the external field. These "loose" spins form 
domains, and between two spins in the same domain there are non-vanishing correlations. 
It was shown then that the connected correlation function %(/) is primerely determined 
by the probability of having a DLS of size / in the system. Finally, this probability was 
calculated using an analogy with the surviving probability of some random walker on a 
finite intervall. 

Considering different types of random field and random bond distributions we have 
calculated the correlation length using the transfer matrix method. Analysing these results 
one my observe two different behaviour in the weak disorder limit. The correlation length 
either vanishes as £ ~ l/log(l/p) or it is divergent like £ ~ 1/p. The former behaviour is 
found in the RFIM with binary distribution in eq(18) as well as for the RBIM in the pure 
ferromagnetic limit q — > in eq(38). On the other hand the correlation length is diverging 
in the diluted RFIM in eqs(20) and (28) as p — > 0. Similar p dependence is observed in 
the RBIM in the antiferromagnetic limit, i.e. when p — > in eq(34). 

The exponential decay of correlations for random Ising chains is found as a general 
rule. One may find, however, a slower decay of correlations if the strength of randomness is 
smoothly position dependent. Let us consider a semi-infinite RFIM in which the strength 
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of the random field decays like hi ~ i~ s from the surface. Then the equivalent random 
walker has to be considered on an intervall the size of which is increasing in time. Using 
results about random walkers [25] and directed polymers [24] inside a parabola, one can say 
that the decay of correlations in this inhomogeneuosly disordered system is of a stretched 
exponential form for < s < 2, whereas it can be described as a power law for s > 2. 
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Figure Caption: 



.1. The integrated random field function H(i) for some fixed values of the random field. 
The corresponding ground state of the system at a given value of the coupling J is 
indicated below. Here vertical dashed lines denote the possible domain wall positions. 
State of spins in the region {7, 11} is not fixed, they form a DLS(see text). 
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The integrated random field function H(i) for some fixed values of the random field. 
The corresponding ground state of the system at a given value of the coupling J is indicated 
below. Here vertical dashed lines denote the possible domain wall positions. State of spins 
in the region {7, 11} is not fixed, they form a DLS(see text). 
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